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I. 


INTRODUCTION 


High  Energy  Laser  (HEL)  technology  has  evolved  to  the  point  that  laser  weapons 
now  have  the  capacity  to  be  used  in  naval  platforms.  Their  advantages  over  conventional 
weapons  in  military  applications  include  unlimited  magazines  (as  long  as  they  are  fed 
with  power),  low  cost  per  shot,  and  very  precise  delivery  of  the  payload  at  the  speed  of 
light. 

However,  unlike  conventional  weapons,  laser  weapon  performance  is  greatly 
affected  by  weather  conditions  and  atmospheric  properties.  Thus,  the  atmospheric 
propagation  of  HELs  needs  to  be  studied  in  order  to  characterize  their  effectiveness  and 
performance  as  laser  weapons. 

In  this  thesis,  the  effects  of  relative  target  and  platform  motion  on  the  atmospheric 
propagation  of  HELs  is  studied,  in  particular  how  this  relative  motion  affects  an 
atmospheric  phenomenon  called  thermal  blooming.  This  study  is  conducted  by  modifying 
a  scaling  code  called  ANCHOR  that  was  developed  by  the  Directed  Energy  Physics 
Group  at  the  Naval  Postgraduate  School  to  include  motions  of  the  platform  and  the  target. 
We  begin  this  study  with  an  overview  of  directed  energy  weapons  in  Chapter  II.  Next,  we 
discuss  the  atmospheric  propagation  of  high  energy  lasers  in  Chapter  III.  In  Chapter  IV, 
we  discuss  how  the  performance  of  the  laser  weapons  is  characterized  and  how  the 
relative  motion  affects  the  atmospheric  propagation.  Then,  we  study  different  cases  in 
order  to  investigate  how  the  propagation  is  affected  by  the  relative  motion  in  Chapter  V. 
Einally,  we  conclude  our  study  Chapter  VI. 
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II.  OVERVIEW  OE  DIRECTED  ENERGY  WEAPONS 


A.  HISTORY 

Directed  energy  (DE)  weapons  are  used  in  order  to  damage  a  target  by  depositing 
energy  onto  it  for  a  specific  dwell  time.  The  most  studied  and  promising  version  of 
directed  energy  weapons  is  lasers  (light  amplification  by  stimulated  emission  of 
radiation). 

The  idea  of  a  laser  started  with  Albert  Einstein  in  1916  as  he  predicted  how 
stimulated  emission  at  particular  wavelengths  can  be  achieved  .  The  first  laser,  the  ruby 
laser,  was  demonstrated  by  Ted  H.  Maiman  in  May  1960,  which  had  around  1  watt  of 
output  power  [1].  After  this  accomplishment,  the  importance  of  this  new  technology  was 
realized  and  further  research  led  to  the  first  gas  laser  (He-Ne  laser)  demonstration  in 
December  1960,  the  demonstration  of  the  first  solid-state  laser  (Nd:YAG)  in  1964,  the 
demonstration  of  the  first  carbon  dioxide  laser  in  1964,  and  the  demonstration  of  the  first 
free  electron  laser  in  1977  [1]. 

The  technology  in  this  area  evolved  to  the  point  that  laser  weapons  now  have  the 
capacity  to  be  used  in  naval  platforms.  A  solid  state  laser  (SSL)  weapon  was  tested  in 
April  2011  on  the  USS  Paul  Foster  (DD  964)  in  the  Pacific  Ocean,  tracking  and  setting 
fire  to  multiple,  small  unmanned  boat  targets.  The  laser  system  was  used  at  high  power 
mode  more  than  35  times  and  withstood  actual  maritime  conditions  including  eight-foot 
waves,  winds  of  25  knots,  and  fog  [2]. 

B.  ADVANTAGES 

In  order  to  understand  the  advantages  of  directed  energy  weapons,  we  need  to 
look  for  the  reason  why  they  are  wanted  to  be  used  in  military  applications  and  what  are 
their  capabilities  that  conventional  weapons  do  not  have. 

The  DE  weapons  apply  damage  at  light  speed  with  very  good  precision.  The 
speed  can  overcome  the  target’s  ability  to  maneuver  kinetic  weapons,  and  the  precision 
enables  damage  to  a  specific  part  of  the  target  that  can  enhance  the  lethality  while  also 
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reducing  the  chance  of  collateral  damage.  The  precision  is  at  the  level  that  the  main  beam 
of  the  Boeing  Airborne  Laser,  a  megawatt-class  chemical  laser,  was  able  to  hit  targets 
that  were  500  kilometers  away  with  pinpoint  accuracy  [3]. 

Another  important  advantage  of  electric  DE  weapons  is  an  essentially  unlimited 
magazine  if  the  platform  can  support  the  weapon  with  power,  unlike  a  conventional 
weapon  that  has  a  specific  amount  of  munition.  This  is  especially  important  in  swarm 
attacks  where  conventional  weapons  may  be  exhausted  before  all  the  targets  are 
eliminated.  Furthermore,  the  cost  per  shot  for  a  DE  weapon  is  cheaper  than  a 
conventional  weapon.  “The  shipboard  fuel  needed  to  generate  the  electricity  for  firing  an 
electrically  powered  laser  would  cost  about  a  dollar  per  shot”  [4].  Finally,  the  effects  of 
the  engagement  on  a  target  can  be  varied  by  changing  the  power  of  the  DE  weapon.  As  a 
result  of  this,  nonlethal  and  lethal  applications  may  be  supported  by  a  single  weapon  [1]. 

C.  TECHNOLOGIES 

1.  Solid  State  Lasers 

An  SSL  uses  a  solid  substrate  that  contains  dopants  as  a  gain  medium.  These 
dopants  are  excited  by  a  pumping  source  in  order  to  achieve  a  lasing  transition.  A  cavity 
that  consists  of  two  mirrors,  one  reflective  and  one  partially  transmittive,  is  used  to 
increase  the  intensity  of  the  light  produced  by  the  lasing  transition  via  the  process  of 
stimulated  emission. 

Two  kinds  of  SSLs  that  will  potentially  be  used  in  military  applications  in  the 
near  term  are  fiber  lasers  and  slab  lasers.  Slab  lasers  have  a  gain  medium  formed  as  a 
rectangular  slab,  while  a  fiber  laser  uses  a  long  cylindrical  medium  consisting  of  an  inner 
and  an  outer  cladding.  Fiber  lasers  have  the  advantage  of  easier  cooling  than  slab  lasers 
due  to  their  geometry,  while  more  light  coupling  can  be  achieved  with  slab  lasers  due  to 
their  large  aspect  ratio. 

Since  they  are  a  more  mature  technology,  and  relatively  compact  compared  with 

Free  Electron  Eases  (FEEs),  SSEs  are  planned  to  be  used  in  DE  weapon  applications  in 

the  near  term.  However,  unlike  FELs,  they  are  restricted  to  operation  at  a  specific 

wavelength,  which  changes  with  the  material  used  as  the  gain  medium.  Also,  due  to  the 
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heat  produced  inside  the  solid  gain  medium  operation,  the  possibility  of  damaging  the 
medium  exists  and  limits  the  output  power  levels  of  SSLs.  FELs  do  not  have  this 
problem.  Also,  in  order  to  reach  the  desired  DE  application  levels  for  SSEs,  multiple 
lasers  need  to  be  combined,  which  results  in  a  decreased  beam  quality  and/or  increased 
system  complexity. 

2.  Free  Electron  Lasers 

Eree  Electron  Easers  (EEEs)  are  used  to  harvest  the  energy  of  unbounded 
relativistic  electrons  as  coherent  light.  There  are  two  main  types  of  EEEs  used:  oscillator 
EEEs  and  amplifier  EEEs.  Either  type  has  several  core  components.  An  injector  and 
accelerator  are  the  first  parts,  where  the  electron  beam  is  created  and  accelerated  to 
relativistic  speed.  Then  these  electrons  pass  through  the  undulator,  where  a  periodically 
changing  magnetic  field  forces  the  electrons  to  wiggle  and  emit  light.  In  an  oscillator 
EEE,  the  light  is  stored  in  a  resonator  and  builds  up  over  many  passes.  In  an  amplifier 
EEE,  the  light  from  a  seed  laser  source  is  amplified  over  a  single  pass.  While  EEEs  are 
inefficient  at  extracting  energy  from  the  relativistic  electrons  on  a  single  pass,  the  leftover 
energy  retained  in  the  electron  beam  can  be  very  efficiently  recovered  to  enhance  the 
overall  wall  plug  efficiency  [1]. 

Since  an  EEE  does  not  have  a  solid  or  gaseous  gain  medium,  thermal  effects  are 
not  as  limiting  a  factor  like  other  types  of  lasers  for  both  beam  quality  and  output  power 
levels.  The  “basic  architecture  of  an  EEE  offers  a  clear  potential  for  scaling  up  to  power 
levels  of  one  or  more  megawatts”  [4].  Oscillator  EEEs  are  not  restricted  to  be  operated  at 
certain  wavelengths,  so  the  wavelength  can  be  tuned  to  match  different  atmospheric 
transmission  sweet  spots.  The  wavelength  of  an  amplifier  EEE  is  fixed  by  the  seed  laser 
used. 

EEEs  are  typically  large,  expensive,  and  the  technology  used  is  not  as  mature  as 
SSEs  [1].  Another  drawback  of  EEEs  is  that  they  create  harmful  radiation  during 
operation  that  needs  to  be  shielded  for  both  personnel  and  material  safety.  This  may  limit 
the  use  of  EEEs  in  military  applications  in  near  term,  but  their  capability  of  reaching  MW 
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power  levels  makes  them  a  viable  choice  for  the  applications  requiring  these  power 
levels. 
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III.  ATMOSPHERIC  PROPAGATION  OF 
HIGH  ENERGY  LASERS 


A.  ATMOSPHERIC  EXTINCTION 

Light  loses  some  of  its  energy  along  the  path  it  transverses  due  to  atmospheric 
extinction,  which  includes  scattering  and  absorption. 

Scattering  occurs  when  photons  interact  elastically  with  particles  in  the 
atmosphere  that  results  in  deflection  of  some  of  the  photons  in  arbitrary  directions.  “The 
particles  responsible  for  scattering  cover  the  sizes  from  gas  molecules  (~10^*cm)  to  large 
raindrops  and  hail  particles  (~1  cm)”  [5].  Absorption  is  an  inelastic  process  in  which 
some  of  the  photons’  energy  is  transferred  to  particles,  and  ultimately  results  in  heating 
the  atmosphere.  This  heating,  in  turn,  can  distort  the  laser  beam  due  to  thermal  blooming. 
Thermal  blooming  will  be  discussed  further  in  Chapter  III,  Section  C. 

The  effects  of  scattering  and  absorption  are  quantified  by  Beer’s  Law, 

P{z)  =  P,e^\  (1) 

In  Equation  (I),  P{z)  is  the  laser  power  after  traveling  a  distance  z  ,  P^  represents  the 

initial  power  at  z  =  0  (i.e.,  living  the  beam  director)  and  s  represents  the  total  extinction 
coefficient  (in  units  of  inverse  length),  which  consists  of  absorption  and  scattering 
coefficients  due  to  molecular  and  aerosol  effects, 

8  =  aa+am+  Pm-  (2) 

In  Equation  (2),  a  and  /?  represent  absorption  and  scattering  coefficients,  while  the  a 
and  m  subscripts  stand  for  aerosol  and  molecular  contributions,  respectively.  Each 
coefficient  in  Equation  (2)  “depends  on  the  wavelength  of  the  laser  radiation”  [6].  In 
order  to  understand  absorption  and  scattering  thoroughly,  aerosol  and  molecular  effects 
need  to  be  discussed  separately. 

I.  Molecular  Effects 

Since  the  wavelengths  used  in  HEEs  are  in  the  infrared  region,  so  all  of  the  effects 
need  to  be  studied  in  this  region  specifically. 
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The  atmosphere  is  composed  of  a  group  of  gases  with  nearly  fixed  concentration 
and  a  group  of  gases  with  variable  concentration  [5].  Nitrogen  (N2)  and  oxygen  (O2) 

compose  ~  99  %  of  the  atmosphere.  However,  they  are  not  responsible  for  the  absorption 
in  near-infrared  region  where  SSLs  operate,  since  they  are  homonuclear  molecules  and 
lack  dipole  moments.  As  seen  in  Figure  1,  water  vapor  ( H2O )  and  carbon  dioxide  ( CO2 ) 

are  predominantly  responsible  for  the  atmospheric  absorption  in  the  infrared  region.  This 
dominance  is  the  result  of  having  many  vibrational  modes.  Molecular  absorption  is 
therefore  a  strong  function  of  the  wavelength. 


0, 


H,0  COt  H|0 

M|0  CO|  0| 

,111  11  11  1 


CO,  0,  H,0  CO,  CO, 

I  Ml 


Figure  1.  Atmospheric  transmittance  measured  over  1820-m  horizontal  path  at 

sea  level.  Source:  [6]. 

When  the  wavelength  of  the  laser  is  considerably  larger  than  the  particle’s  sizes, 
the  scattering  process  is  described  by  Rayleigh  scattering.  The  Rayleigh  scattering  cross 
section,  which  is  proportional  to  the  likelihood  of  a  scattering  event  occurring,  can  be 
approximated  as  [5] 
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2 


(3) 


where  a  is  the  diameter  of  the  molecules  and  assumed  to  be  much  smaller  than  the 
wavelength  of  the  laser  beam  (a  <sc  2) ,  and  n  is  the  refractive  index  of  the  molecules.  As 
wavelength  increases,  the  effect  of  Rayleigh  scattering  on  atmospheric  propagation 
decreases  rapidly  (cr^  oc  .  The  sky’s  blue  color  during  daytime  can  be  explained  by 

Rayleigh  scattering;  since  blue  light’s  wavelength  is  shorter  than  red  light’s,  blue  light  is 
scattered  more  than  red  light. 


Using  the  Rayleigh  scattering  cross  section,  the  molecular  scattering  coefficient 
can  be  calculated  as  [5] 

/^JZ)  =  aM)]N(r)dr,  (4) 

where  A(r)  is  the  number  density  of  molecules  with  radius  r . 


2.  Aerosol  Effects 

Aerosol  extinction  coefficients  are  less  dependent  on  the  light  wavelength  than 
molecular  extinction  coefficients.  These  coefficients  are  determined  by  the  combined 
effects  of  the  aerosol’s  index  of  refraction  and  the  size  distribution  of  aerosol  particles, 
according  to  the  theory  of  Mie  scattering  [7].  In  general,  the  aerosol  effects  are  more 
prominent  at  shorter  wavelengths,  at  least  in  the  IR  spectral  region  [7]. 

The  effect  of  aerosol  scattering  on  atmospheric  propagation  is  generally  far 
greater  than  molecular  (Rayleigh)  scattering  except  at  very  short  wavelength  values  due 
to  molecular  scattering’s  dependence  on  wavelength  as  2 Figure  2  shows  that  for 
wavelengths  up  to  ~  0.3  pm ,  molecular  (Rayleigh)  scattering  dominates  the  total 
scattering  extinction,  but  aerosol  scattering  starts  to  become  dominant  at  near-IR  and 
longer  wavelengths. 
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Figure  2.  Calculated  aerosol  absorption  and  extinction  coefficients  for  23-km 
visibility  conditions  and  a  continental  aerosol  model.  Source:  [8]. 


B.  ATMOSPHERIC  TURBULENCE 

Atmospheric  turbulence  is  the  result  of  temperature  differences,  convection,  wind 
shear,  and  inertial  cascades  of  larger  scale  turbulences  to  smaller  scales  [1]. 
Inhomogeneities  in  the  atmosphere  (eddies)  result  in  density  variations  and,  in  turn,  index 
of  refraction  differences  along  the  laser  beam  path.  “The  shimmering  of  distant  objects  on 
a  hot  day  or  the  twinkling  of  the  stars  at  night”  are  examples  caused  by  index  of 
refraction  variances  due  to  atmospheric  turbulence  [9].  These  index  differences  also 
distort  the  wave  fronts  of  an  HEL,  potentially  causing  the  laser  beam  to  break  up  into 
beamlets  that  travel  in  slightly  different  directions. 

An  important  parameter  that  is  used  to  characterize  the  severity  of  atmospheric 
turbulence  is  the  refractive  index  structure  coefficient  with  units  of  m^^^^ .  This 
parameter  increases  as  the  temperature  differences  of  the  eddies  increases,  so  an  increase 
in  cl  means  an  increase  in  atmospheric  turbulence.  Cl  tends  to  decrease  with  altitude. 

Typical  Cl  values  near  the  surface  vary  from  about  10 for  relatively  weak 
turbulence  to  10  m  for  relatively  strong  turbulence.  Furthermore,  turbulence  varies 
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throughout  the  day  and  tends  to  be  smallest  “about  one  or  two  hours  before  sunrise  and 
after  sunset,”  while  the  peak  values  tend  to  be  measured  around  noon  on  sunny  days  [9] . 
The  irradiance  patterns  for  a  laser  propagating  for  5  km  through  a  region  of  constant 
are  shown  in  Figures  3  through  5.  In  Figure  3,  the  turbulence  is  extremely  weak  at 
=  10 ,  which  has  almost  no  effect  on  the  irradiance. 


X  dimension  [m] 

2  _ 10  _ 2/3 

Figure  3.  Irradiance  with  weak  turbulence  (C„  =10  m  '  ).  Source:  [10]. 

The  irradiance  pattern  with  stronger  turbulence  (C^  =10 is  shown  in 

Figure  4.  The  figure  shows  that  the  laser  beam  gets  distorted  as  the  atmospheric 
turbulence  increases. 
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Figure  4.  Irradiance  with  moderate  turbulenee  ( =  10 ). 

Source:  [10]. 


The  irradiance  pattern  with  even  stronger  turbulence  (C^  =5x10  m  is 
shown  in  Figure  5.  The  figure  shows  that  the  distortion  level  reaches  to  the  point  where 
the  peak  irradiance  is  down  by  over  an  order  of  magnitude  relative  to  Figure  3. 
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Figure  5.  Fradiance  with  strong  turbulence  ( =  5  x  10 m^^^^ ). 

Source:  [10]. 
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The  Fried  parameter  is  the  diameter  over  which  transverse  beam  coherence  is 
maintained.  It  incorporates  information  about  along  the  path,  the  wavelength  of  the 
laser  and  the  propagation  distance.  Larger  values  of  correspond  to  weaker  turbulence 
and  more  complete  constructive  interference  at  the  target.  For  a  focused  beam 
propagating  along  a  path  of  constant  ,  the  Fried  parameter  is  given  by  [1] 


Tq  =  0.33 


/I'’/' 


3/5 


(5) 


where  A  is  the  wavelength  of  the  laser  light  and  R  is  the  range  of  the  target.  Larger 
values  and  longer  ranges  result  in  smaller  revalues,  while  longer  wavelengths  result  in 
longer  revalues.  Generally,  turbulence  becomes  significant  when  is  smaller  than  the 
beam  director  size.  As  an  example,  the  Fried  parameter  is  « 12.5  cm  for  A  =  l  pm, 
cl  =  10 and  7?  =  5  km . 


C.  THERMAL  BLOOMING 

Absorption  of  laser  light  by  the  molecules  and  aerosols  in  the  atmosphere  results 
in  heating  the  air,  which  then  changes  the  density  of  air  along  the  laser  path.  These 
density  gradients  create  refractive  index  changes  in  this  path.  These  alterations  in  the 
optical  properties  of  the  path  cause  the  air  to  act  like  a  concave  lens,  which  in  turn  results 
in  defocusing  and  spreading  of  the  laser  beam.  This  nonlinear  atmospheric  phenomena  is 
called  thermal  blooming.  The  effect  of  thermal  blooming  is  especially  crucial  in  low 
altitude  engagements  (due  to  higher  water  vapor  and  aerosol  content)  such  as  ground/sea 
surface-based  scenarios,  as  well  as  slow  moving  targets  in  light  winds  due  to  reduced 
beam  slewing.  Wind  across  the  beam  path  tends  to  reduce  the  severity  of  thermal 
blooming  since  it  replaces  the  air  warmed  by  the  laser  with  cooler  air. 

There  are  two  limiting  cases  of  thermal  blooming:  conduction  dominated  thermal 
blooming  and  convection  dominated  thermal  blooming  [11].  Conduction-dominated 
thermal  blooming  exists  for  the  cases  where  there  is  basically  no  wind  (v=0).  This  case 
is  rarely  met  for  the  HEL  applications  along  the  entire  beam  path,  though  there  can  be 
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places  along  the  path  that  are  stagnant.  Convection-dominated  thermal  blooming  exists 
when  there  is  the  wind  or  beam  motion.  This  case  is  met  for  most  HEL  applications  and 
has  been  studied  most  extensively  [11]. 

When  wind/beam  motion  exists,  cooler  air  enters  from  one  side  of  the  beam  and 
warms  as  it  transverses  the  beam.  This  causes  the  air  to  heat  up  in  the  downwind 
direction,  and  thus  the  refractive  index  drops  in  the  downwind  direction.  Since  light 
bends  toward  higher  refractive  indices,  the  irradiance  bends  into  the  wind,  often  forming 
a  distinct  crescent  shape  as  illustrated  in  Figure  6. 


Figure  6.  Typical  thermal  blooming  effected  irradiance  pattern  with  wind 
coming  from  left  side  {P  =  100  kW  ,  R  =  5  km ,  =  10  m / s )  . 

Source:  [10]. 
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IV.  METHODOLOGY 


ANCHOR  (Atmospheric  NPS  Code  for  HEL  Optical  pRopagation),  a  scaling 
code  developed  by  the  NPS  DE  group,  is  used  to  characterize  laser  performance.  The 
code  incorporates  parameters  like  turbulence,  wind,  jitter,  and  extinction.  Eor  this  thesis, 
the  effect  of  relative  target  and  platform  motion  on  atmospheric  propagation  was  also 
added  to  the  code.  Eor  different  laser  power  and  laser  wavelength  inputs,  time  averaged 
irradiance  and  power  in  the  bucket  plots  can  be  acquired  using  this  code.  The  code  takes 
atmospheric  parameters  from  other  models  like  EEEDR,  NAVSEAM,  MODTRAN,  and 
CO  AMPS. 

In  order  to  estimate  the  irradiance  at  the  target,  a  couple  of  parameters  need  to  be 
taken  into  account.  Atmospheric  extinction,  beam  quality,  platform  jitter,  atmospheric 
turbulence,  and  thermal  blooming  are  all  effects  that  cause  the  irradiance  to  decrease 
along  the  propagation  path.  Assuming  a  Gaussian  beam  at  the  source,  the  time-averaged 
irradiance  profile,  with  platform  jitter  and  turbulence-induced  wander  and  spreading 
effects  included,  is  estimated  as  [7] 

(6) 

where  is  the  total  power  leaves  the  beam  director,  s  is  the  total  atmospheric 
extinction  coefficient,  I  is  the  range  to  the  target,  is  the  time-averaged  1/e  radius  of 
the  mean  irradiance  profile  in  the  absence  of  thermal  blooming,  Sj^  represents  the  Strehl 
ratio  that  accounts  for  the  effect  of  the  thermal  blooming  on  the  time  averaged  irradiance. 
The  time-averaged  beam  radius  on  the  target  is  given  by 


In  Equation  (7),  the  j ,  d  and  t  subscripts  represent  contributions  from  beam  jitter, 
diffraction  and  atmospheric  turbulence,  respectively. 
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A.  LINEAR  EFFECTS  ON  SPOT  SIZE 
I.  Beam  Jitter 

Beam  jitter  effect  is  caused  by  vibrations  and/or  tracking  errors  on  the  beam 
director  platform.  The  time  averaged  mean  square  radial  displacement  of  the  focal  spot 
due  to  beam  jitter  in  Equation  (7)  is  given  approximately  as  [7] 

~  (8) 
where  6^^  represents  the  angular  variance  due  to  jitter  and  i  is  the  range  to  the  target. 
Typical  values  of  angular  jitter  are  on  the  order  of  a  few  jurad  . 


2. 


Diffraction 


The  effect  of  diffraction  on  the  laser  beam  spot  size  for  an  infinite  Gaussian 
source  can  be  estimated  as  [7] 


k  Wn 


+  w„ 


A 


1-  — 

V  Fj 


(9) 


where  i  represents  the  range  to  the  target,  k  =  Itt/X  is  the  wave  number,  is  the  He 

beam  radius  at  the  beam  director,  F  is  the  focal  range  of  the  laser  beam,  and 
represents  the  beam  quality  factor.  The  value  is  the  measurement  of  how  many  times 
an  ideal  beam  can  fit  across  the  real  beam  in  each  transverse  direction  at  the  focus.  For  an 
ideal  Gaussian  beam,  M^=l  whereas  >1  for  non-ideal  laser  beams  [12].  For  this 
thesis,  the  focal  range  F  is  equal  to  the  target  range  £ .  For  this  condition,  the  second 
term  in  the  square  root  drops  out  and  Equation  (9)  simplifies  to 

M£ 


= 


kw„ 


(10) 


3.  Turbulence 

The  effect  of  the  atmospheric  turbulence  (see  Chapter  III,  Section  B)  on  the  laser 
beam  spot  size  can  be  estimated  as  [13] 


11 

w,  «— ,for  r^<D, 
kr^ 


(11) 
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where  is  the  Fried  parameter  (see  Equation  (5))  and  D  represents  the  diameter  of  the 
beam  director.  The  effect  of  atmospheric  turbulence  increases  with  the  range  of  the  target 
and  decreases  with  the  Fried  parameter  since  it  represents  the  diameter  that  the 
transverse  coherence  of  the  laser  beam  is  maintained. 


B.  THE  EFFECT  OF  THERMAL  BLOOMING  ON  SPOT  SIZE 


The  effect  of  the  thermal  blooming  on  the  laser  beam  propagation  through  the 
atmosphere  is  discussed  in  Chapter  III,  Section  3.  The  effect  on  the  laser  beam  spot  size 
can  be  estimated  mathematically  using  the  thermal  blooming  distortion  number,  which  is 
calculated  as  [1] 

fMilIhKhffc,  (12) 

p  path  ^wind  (z)D(z)  ’ 

where  z  is  the  distance  along  the  beam  path,  P  is  the  initial  power  of  the  laser,  is  the 


ambient  air  density,  is  the  specific  heat  at  constant  pressure,  a  is  the  absorption 
coefficient,  T  is  the  transmission  of  the  laser  energy,  n^.  =dnldT  =-(n()-l)/r  is  the 
change  in  the  refractive  index  of  air  with  respect  to  temperature  where  italic  T  represents 
the  temperature,  represents  the  effective  wind  speed  perpendicular  to  the  laser  beam 


and  D  is  the  beam  diameter.  Absorption  near  the  beam  director  causes  more  significant 
thermal  blooming  distortion  than  the  absorption  near  the  target  [1].  This  effect  is  called 
lever  arm  effect.  In  order  to  compensate  this  effect,  the  following  weighting  function  can 
be  included  in  the  integral  in  Equation  (12)  [1], 


F{z)  =  - 


1  + 


T25^ 
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z/i 
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where  D  represents  the  primary  aperture  diameter  in  this  case. 


(13) 


If  there  is  no  change  in  the  atmospheric  conditions  along  the  path.  Equation  (12) 
can  be  simplified  as  [1] 

-AyflkPriT-aT  . 
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(14) 


Thermal  blooming  effects  become  significant  when  >25  [1];  the  laser  beam  gets 
more  distorted  as  increases.  In  Equation  (14),  it  is  seen  that  the  transverse  effective 
wind  decreases  the  distortion  level  due  to  thermal  blooming,  and  in  turn  can  increase  the 
irradiance  at  the  target.  The  irradiance  patterns  for  a  laser  (P  =  100  kW  )  propagated  5 
km  through  a  region  with  different  effective  wind  speeds  are  shown  in  Figures  7  and  8.  In 
Figure  7,  the  effect  of  the  thermal  blooming  on  irradiance  is  shown  for  =3  m  /  s  . 

Irradiance  xio* 

12 

10 

8 

CM 

E 

6  ^ 

4 

2 

0 

-0.4  -0.2  0  0.2  0.4 

X  dimension  [m] 

Figure  7.  Thermal  blooming  effected  irradiance  patterns  with  =  3  m  /  s 

=465,P  =  100  kW,  P  =  5  km).  Source:  [10]. 

In  Figure  8,  the  effect  of  the  thermal  blooming  on  irradiance  is  shown  for 
'^wind  =5  m/s .  It  is  seen  that  as  the  transverse  wind  speed  increases,  the  irradiance  also 

increases  (note  the  different  irradiance  scales  on  the  graphs),  although  the  crescent¬ 
shaped  distortion  of  the  beam  is  enhanced. 
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Figure  8.  Thermal  blooming  effected  irradiance  patterns  with  =  5  m  /  s 
{N^=219  ,P  =  \00  kW,  i?  =  5  km).  Source:  [10]. 


The  Strehl  ratio  Sj.^,  the  ratio  of  the  peak  irradiance  of  the  thermal  blooming 
distorted  beam  to  the  peak  irradiance  of  a  reference  beam  with  no  thermal  blooming  can 
be  calculated  as  [1] 


5  ^ 

™  h  l  + 

where  a  =  0.0625  and  m  =  2  for  a  Gaussian  distributed  beam,  whereas  a 
m  =  1.2  for  a  uniform  distributed  laser  beam. 


(15) 

0.01  and 


Thermal  blooming  also  affects  HEL  laser  systems  such  that  after  a  certain  laser 
output  power,  ,  the  irradiance  on  the  target  starts  to  decrease.  This  critical  power 

level  can  vary  with  the  laser  wavelength  and  the  atmospheric  and  weather  conditions,  and 
can  be  estimated  by  [1] 


P  = 

crit 


f  p  ^ 

r  -1  / 

UJ 

1/m 


(16) 
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In  ANCHOR,  the  distortion  number  from  Equation  (14)  is  calculated  as  well  as  a 
weighted  distortion  number  which  uses  the  same  equation  but  includes  a  weighting 
function  granting  more  effect  of  thermal  blooming  to  the  midpoint  along  the  beam  path. 
Using  these  two  numbers,  an  analytic  expression  for  the  Strehl  ratio  is  used  that  is 
calibrated  against  full  diffraction  codes. 

1.  Effective  Transverse  Wind  Speed 

In  order  to  calculate  the  effective  transverse  wind  speed  across  the  beam,  both 
target  and  platform  motions  and  the  natural  wind  need  to  be  taken  into  consideration.  As 
a  result  of  the  relative  target  and  platform  motion,  the  beam  director  must  slew  to  keep 
the  laser  on  the  target.  This  slew  motion  also  introduces  an  effective  wind  across  the 
beam. 

In  Figure  9,  the  parameters  needed  to  calculate  the  transverse  wind  speed  is  given. 
The  calculations  are  made  relative  to  the  ship  using  spherical  coordinates. 


The  target  velocity  relative  to  the  ship  is  calculated  as 


V  =V  -V, 

target  target  ship  ’ 


(17) 


where  represents  the  target  velocity  and  is  the  velocity  of  the  ship.  Likewise, 


the  wind  velocity  relative  to  the  ship  is  calculated  as 
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(18) 


^wind  ^wind  ^ ship  ’ 


where  represents  the  natural  wind  velocity.  Angular  velocity  of  the  target  about  the 
ship  is  calculated  using  the  relative  target  velocity  as 


CO  = 


r  xV 

target  target 


r 

target 


(19) 


where  represents  the  displacement  vector  of  the  target  relative  to  the  ship.  Slew 
velocity  along  path  r  is  calculated  using  the  angular  velocity  of  the  target  as 


Vsie.=^xr. 


(20) 


The  effective  wind  velocity  is  calculated  using  the  slew  velocity  and  relative  wind 
velocity  as 

yeJf=Kind-ysle.-  (21) 

Finally,  the  effective  transverse  wind  (i.e.,  the  cross  wind),  is  calculated  using  the 
effective  velocity  as 


=(%4)h{%-e)0. 


where  (j)  and  6  are  the  spherical  unit  vectors,  so 

=  -  sin  (l)x  +  co^(l)% 
and 


(22) 


(23) 


0  =  cos  0 cos  x+  cos  0  sin  (jdy  -  sin  Oz.  (24) 

We  created  a  MATLAB  script  using  these  equations  that  plots  the  effective 
transverse  wind  speed  for  a  constant  target  range  and  varying  orientation  about  a  circle 
centered  on  the  ship  with  the  same  and  parameters.  In  Figure  10,  the 

effective  transverse  wind  speed  is  plotted  for  =10(m/s)y  and 


ywind  =-10(m/s)y  .  In  these  figures,  arrow  length  represents  the  magnitude  of  the  cross 

wind  along  the  beam.  Also,  the  color  of  the  arrows  evolves  from  blue  to  yellow  to  red  as 
the  magnitude  increases.  In  this  example,  there  is  no  cross  wind  along  the  vertical  axis, 
because  in  that  case  the  wind  is  in  the  same  direction  as  the  ship  and  target  motion,  so  a 
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stagnation  zone  exists  in  this  direction.  This  region  is  where  thermal  blooming  will  be  the 
most  problematic.  Along  the  horizontal  axis,  the  effect  of  the  thermal  blooming  on 
propagation  is  the  least  because  the  strongest  cross  wind  exists  along  this  path. 


Figure  10.  Effective  transverse  wind  speed  for  =  10(m  /  s)  y  and 

^v.,w=-10(m/s)^ 


In  Figure  1 1,  the  effective  transverse  wind  speed  is  plotted  for  =  10(m  /  s)  y , 

^target  =  10(tTi  /  s)  X  uud  =  10(m  /  s)x  +  10(m  /  s)y .  Here,  we  also  see  that  a  stagnation 

zone  exists  for  certain  target/ship  orientations.  Furthermore,  the  position  along  the  path 
where  this  stagnation  zone  occurs  changes  with  the  orientation. 
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Figure  1 1 .  Effective  transverse  wind  speed  for  =  10(m  /  s)  y , 
V,arge,  =  /  s)  X  and  =  10(m  /  s)x  +  10(m  /  s)y . 


C.  OTHER  PARAMETERS  NEEDED  TO  CHARACTERIZE  LASER 

PERFORMANCE 

Other  than  the  irradiance,  laser  performance  can  be  characterized  by  power  in  the 
bucket  and  dwell  time. 


1.  Power-in-the-Bucket 


The  power  delivered  to  a  circle,  or  “bucket,”  of  specified  radius  on  the  target  is 
called  power-in-the-bucket.  The  calculation  of  the  power  in  the  bucket  changes  with  the 
chosen  bucket  radius  r^in  comparison  with  the  thermal  blooming  affected  beam  radius. 

The  thermal  blooming  affected  beam  radius  can  be  calculated  using  the  Strehl  ratio  as 
[14] 


w, 


fot.eff 


'TB 


(25) 


Since  ^  accounts  for  the  thermal  blooming  effect,  the  irradiance  on  the  target  can  be 


rewritten  as  [14] 


23 


(26) 


(I)= - 

If  ,  power  in  the  bucket  can  be  calculated  approximately  as  [14] 

(27) 

as  all  of  the  delivered  power  falls  within  the  defined  bucket.  If  ,  then  the  power 

in  the  bucket  can  be  calculated  approximately  as  [14] 

PB^{I)[7rr^).  (28) 

In  this  case,  we  are  not  interested  in  all  the  power  delivered  to  the  target.  The  total 
power  delivered  to  the  target  is  more  than  the  power  in  the  bucket. 

2.  Dwell  Time 

A  laser  weapon  is  used  to  transfer  energy  to  the  target,  and  in  turn  damage  it  by 
melting  the  material  used  on  the  target.  The  time  needed  for  melting  through  the  target  is 
called  dwell  time. 


a.  Power  Loss  Mechanisms 


In  order  to  calculate  the  dwell  time,  power  loss  on  the  target  needs  to  be  known. 
There  are  two  power  loss  mechanisms:  conductive  and  radiative  losses.  These  losses 
decrease  the  power  level  transferred  to  the  target. 


Conduction  loss  exists  as  a  result  of  the  heat  increase  at  a  specific  area  on  the 
target  and  the  conduction  of  the  heat  through  the  other  parts  of  the  target.  Conduction  loss 
can  be  quantified  using  Fourier’s  heat  conduction  law  as 


P  = 

cond 


I^Kond  -Tq) 


Ax 


(29) 


where  k  is  the  thermal  conductivity  of  the  material,  is  the  melting  temperature  of  the 


material,  is  the  ambient  temperature  of  the  material,  is  the  cross-sectional 

surface  area  of  the  surrounding  material  and  Ax  is  the  distance  of  the  temperature 
gradient  in  material.  This  distance  depends  on  the  material  properties  and  the  dwell  time; 
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for  typical  metals  and  dwell  times  of  r^„^;/~10s,  A;c~2cm  [15].  Typically,  the 
conductive  power  loss  is  much  greater  than  the  radiative  power  loss. 

Radiative  power  loss  is  due  to  blackbody  radiation  as  the  material  heats  up  and 
glows.  Using  Stefan-Boltzmann  law,  the  power  radiated  away  is  proportional  to  the 
fourth  power  of  the  blackbody’ s  temperature  azT^)  ,  the  radiated  power  due  to 
temperature  change  after  laser  dwell  can  be  calculated  as 

=  Aea{T: -T:),  (30) 

where  A  is  the  area  over  which  the  radiation  occurs,  ^is  the  material  emissivity  and  cr 
is  the  Stefan-Boltzmann  constant  ( ~  5.670 x  10  *  W  m'^  K'"* ). 


Total  power  loss  due  to  conduction  and  radiation  is  calculated  as 


p  =  P  +  P 

loss  cond  rad. 


(31) 


The  power  available  to  melt  the  target  can  be  quantified  using  the  power  loss  as 

p  —  p  X  f  —  P  1321 

melt  trans  J  target  ^  loss,  ^  * 


where  represents  the  power  transferred  to  the  target  and  is  the  fractional 

target  absorption.  For  aluminum,  is  ~  0.2.  In  order  to  melt  the  target,  the  power 

transferred  to  the  target  must  exceed  the  power  loss  on  the  target,  otherwise  the  target  can 
not  be  melted. 


b.  Energy  Needed  to  Melt  the  Target 

The  energy  needed  to  melt  a  target  varies  with  the  melting  temperature  of  the 
target  material,  heat  of  fusion  of  the  target  material,  and  the  mass  of  the  material.  The 
energy  needed  to  reach  the  melting  point  of  the  material  is  calculated  as 

=  CpinAT,  (33) 

where  Cp  represents  the  specific  heat  capacity  of  the  material,  m  is  the  mass  of  the 
material  needed  to  melt  and  AT  =T^-T^^  is  the  change  in  the  temperature.  The  energy 
needed  to  melt  the  material  at  the  melting  point  is  calculated  as 

Q^=mNH,  (34) 
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where  A//  represents  the  heat  of  fusion  of  the  target  material.  The  total  energy  to  melt 
the  target  material  is  calculated  as 

Qn,eU=Q\+Q2-  (35) 


c.  Dwell  Time  Calculation 

Dwell  time  can  be  calculated  using  these  calculations  as 

_  _  Qmeh 

dwell  p 

melt 


(36) 


Dwell  time  can  vary  with  the  power  level  transferred  to  the  target,  the  power  loss  on  the 
target,  material  of  the  target,  and  the  volume  needed  to  melt. 


d.  Laser  Performance  Calculation  Example 

As  an  example,  required  on-target  irradiances  for  melting  a  100  cm^  circular  area 
on  an  aluminum  target  are  calculated  for  dwell  times  of  3  and  6  seconds  and  material 
thicknesses  of  1  and  3  mm.  Table  1  shows  the  physical  properties  of  aluminum  used  in 
these  calculations. 


Table  1.  Physical  Properties  of  Aluminum 


Parameter 

Value 

Density 

2.70 

cm 

Thermal  conductivity,  k 

W 

237 

mK 

Distance  of  the  temperature  gradient,  Ax 

~  2  cm 

Melting  temperature, 

933.47  K 

Material  emissivity,  s 

0.05 

Specific  heat  capacity  of  the  material,  Cp 

897  ^ 

kgK 

Heat  of  fusion  of  the  material.  Mi 

400  H 

kg 
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Calculated  total  energies  to  melt  the  target  and  total  power  losses  for  1  mm  and  3 
mm  thicknesses  are  shown  in  Table  2. 


Table  2.  Total  Power  Loss  and  Total  Energy  to  Melt  Calculation  Results 


Parameter 

For  1  mm 
thickness 

For  3  mm 
thickness 

Total  energy  to  melt  the  target, 

26.14kJ 

78.42  kJ 

Total  power  loss, 

2.68  kW 

8.04  kJ 

Using  the  parameters  in  Table  2,  the  required  on-target  irradiances  are  shown  in 
Table  3. 


Table  3.  Required  On-Target  Irradiances  to  Melt  the 
Target  Calculation  Results 


Parameter 

For  1  mm 
thickness 

For  3  mm 
thickness 

Required  on-target  irradiance  to  melt 
target  for  =  3  seconds 

MW 

1.13  , 

m 

3,4  T 

m 

Required  on-target  irradiance  to  melt 
target  for  =  6  seconds 

0.7  “'f 
m 

.  ,  MW 

2.1  — 
m 
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V.  RESULTS 


A.  NOTIONAL  LASER/ATMOSPHERIC  PARAMETERS 

In  order  to  evaluate  the  effects  of  relative  target  and  platform  motion  on  the 
atmospheric  propagation  of  HELs,  the  parameters  used  as  inputs  to  ANCHOR  are  given 
in  Table  4. 


Table  4.  The  Input  Parameters  Used  in  ANCHOR. 


Laser  Output  Power  P 

200  -  1000  kW 

Laser  Wavelength  X 

1.064  pm 

Beam  Quality 

3 

Platform  Height  Hp 

10  m  Above  Sea  Level 

Target  Height  Hp 

0.4  -  1300  m  Above  Sea  Level 

Diameter  of  Beam  at  Source  D 

0.3  m 

Beam  Shape  at  Source 

Uniform 

Beam  Jitter  6^^ 

5x10 rad 

Speed  of  Target 

0  ^  198  m/s 

Speed  of  Platform  1/ 

0^10  m/s 

Speed  of  Wind 

0^3  m/s 

Bucket  Radius 

0.05  m 

Fractional  Target  Absorption 

0.15 
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In  order  to  evaluate  the  extinction  coefficients  for  X  =  1.064  |jm  along  the  laser 
beam  path,  MODTRAN  is  used.  MODTRAN  is  a  commonly  used  radiative  transfer 
model  developed  by  Spectral  Sciences  Inc.  and  the  U.S.  Air  Force  Research  Laboratory 
[16].  The  parameters  used  as  input  for  MODTRAN  are  given  in  Table  5. 


Table  5.  Tropical  Maritime  Input  Parameters  Used  for  MODTRAN. 


Parameter 

MODTRAN  Model 

Geographical-seasonal 
Model  Atmosphere 

Tropical  atmosphere  (15  ° 
Latitude) 

Aerosol  Model 

Maritime  extinction 

Cloud  and  Rain  Model 

No  rain 

Meteorological  Range 

23  km 

Temperature, 

300  K  at  sea  level  (varies  with 
altitude) 

Absolute  Humidity 

19.7  g/m^  (80%  relative  humidity) 
at  sea  level  (varies  with  sea  level) 

The  refractive  index  structure  coefficient  along  the  laser  beam  path  is 
calculated  using  the  Hufnagel- Valley  model.  This  model  characterizes  the  variations  of 
cl  with  the  altitude  for  a  given  surface  value  of  C,J  and  is  specified  by  the  following:  [1] 


C'(/z)  =  5.94x10 


-53 


IT 

V27y 


10  -A/IOOO 


1  lu 

n  e 


-16  -^1500 


+  Ae 


-/i/lOO 


(37) 


where  h  is  the  height  in  meters,  W  is  the  high  altitude  wind  speed  in  m/s  and  A  is  the 
cl  surface  value.  In  this  thesis,  A  is  set  to  the  moderate  value  of  10 m^^^^ . 
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B.  COMPARISON  WITH  WAVETRAIN 


1.  Wavetrain 

In  order  to  validate  the  ANCHOR  code  results,  Strehl  ratio  outputs  are  compared 
with  the  output  of  WaveTrain  for  the  same  conditions.  WaveTrain  is  a  software  tool 
developed  by  MZA  Associates  Corporation  “for  high  fidelity  modeling  of  advanced 
optical  systems  such  as  laser  weapons  systems”  [17].  It  is  a  full  diffraction  code,  which 
potentially  makes  it  more  accurate  than  ANCHOR,  but  ANCHOR  is  much  faster  than 
WaveTrain.  All  of  the  WaveTrain  outputs  for  this  thesis  were  provided  by  Dr.  Conor 
Pogue. 

In  WaveTrain,  the  Strehl  ratios  are  calculated  using  the  same  absorption 
coefficients,  target,  platform  and  wind  velocities,  target  and  platform  heights,  and  target 
distances  with  ANCHOR.  Strehl  ratio  values  are  the  division  of  the  peak  irradiance  with 
thermal  blooming  and  the  peak  irradiance  without  thermal  blooming  for  WaveTrain. 

2.  Results  and  Discussion 

In  Figures  12  and  13,  the  comparison  of  Strehl  ratios  of  ANCHOR  and  Wavetrain 
is  plotted  in  polar  coordinates  for  a  constant  target  range  and  varying  orientation  with  the 

same  and  parameters.  Atmospheric  turbulence  effects  are  ignored  for 

both  codes.  In  Figure  12,  the  parameters  used  are  =  10(m/  s)  y ,  =  10(m/  s)  y , 

= -10(m/s)y  ,  /’  =  100kW, /l  =  1.064  pm  (typical  wavelength  used  in  SSLs), 
R  =  5  km  and  height  of  the  target  hj.  =  10  m.  The  effective  transverse  wind  speed  plot  for 

the  same  platform,  target  and  wind  velocities  was  shown  in  Figure  10.  It  is  seen  that, 
ANCHOR  tends  to  underestimate  the  Strehl  ratio  more  as  the  cross  wind  along  the  laser 
beam  increases. 
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Figure  12.  Strehl  ratio  output  eomparison  of  Wavetrain  and  ANCHOR  for 
A,  =  10(m  /  s)  y ,  V',,,,,  =  10(m  /  s)  y ,  V.,„  =  -10(m  /  s)y  ,  />  =  100 
kW,  X  =  1.064  pm ,  R  =  5  km  and  hj  =  10  m.  Source:  [10]. 

In  Figure  13,  the  parameters  used  are  =10(m/s)y ,  =10(m/s)x, 

=10(m/s).r  +  10(m/s)5',  P  =  100kW,  A  =1.064  pm,  R  =  5km  and  height  of  the 
target  hj,  =  10  m.  The  effective  transverse  wind  speed  plot  for  the  same  platform,  target 

and  wind  velocities  was  shown  in  Figure  1 1 .  It  is  seen  that  the  Strehl  ratios  of  the  two 
methods  agree  generally  except  the  regions  where  ANCHOR  slightly  underestimates  the 
Strehl  ratio  in  the  case  where  the  cross  wind  speed  is  larger  at  the  zone  closer  to  the 
platform  and  slightly  overestimates  the  Strehl  ratio  where  the  cross  wind  is  larger  at  the 
zone  closer  to  the  target. 
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Figure  13.  Strehl  ratio  output  comparison  of  Wavetrain  and  ANCHOR  for 
y,up  =  10(m /  s)  y ,  =  10(m / s) x ,  =  10(m /  s)x  +  10(m /  s)y , 

P  =  100  kW,  Z  =  1.064  pm ,  R  =  5  km  and  =  10  m.  Source:  [10]. 


C.  ANCHOR  EXAMPLES 

The  irradiance  of  the  laser  will  be  plotted  for  a  constant  target  altitude,  hj ,  for 
varying  target  ranges  from  the  platform  at  the  origin.  Every  point  except  the  origin 
corresponds  to  a  target  with  the  same  velocity.  The  velocity  of  the  wind  also 
does  not  change  for  any  target  orientation.  In  Figure  14,  these  target  orientations 

are  displayed  for  several  points  on  the  polar  plot  with  the  velocities  of  =10(m/s)5', 
y,ar,e,  =  10(m /  s) X  and  =  10(m  / s) x  +  10(m / s)  y . 
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Figure  14.  Visualization  of  target  orientation  for  =  10(m  /  s)  y, 
V.arser  =  10 ( iTi / s ) X  and  =10(m/s)x  +  10(m/s)3). 


1.  Simple  Cases 

In  all  of  the  irradiance  plots,  latitudinal  range  is  along  the  y-axis,  whereas 
longitudinal  range  is  along  the  x-axis.  The  contour  lines  on  the  plot  connect  points  of 
constant  irradiance  (outer  contour  line  is  for  1  MW/m^  ,  middle  contour  line  is  for 
5  MW/m^  and  inner  contour  line  is  for  10  MW/m^ )  .  The  magnitude  of  the  on-target 
irradiance  is  defined  with  the  color  bar  scale.  In  Figure  15,  the  on-target  irradiance  is 
plotted  for  =  -3(m  / s)y,  hj  =10m  and  P  =  200 kW.  There  is  no  relative  platform 

and  target  motion  for  this  plot.  Since  we  have  a  stagnation  zone  on  the  vertical  axis,  and 
the  cross  wind  gets  stronger  along  the  horizontal  axis,  we  have  an  on-target  irradiance 
pattern  that  is  compressed  along  the  vertical  axis  and  elongated  along  the  horizontal  axis. 
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-5  0  5 

Longitudinal  Range  (km) 

Figure  15.  Simple  case  scenario  for  =  -3(m  /  s)  hj  =10  m,  P  =  200  kW 

and  /I  =  1.064  |jm. 

In  Figure  16,  the  on-target  irradiance  is  plotted  for  =  10(m/ s) 
Kin,  =-3(m  /  s)y,  hj  =  10  m  and  P  =  200  kW.  With  the  introduction  of  platform  motion, 

the  magnitude  of  the  cross  wind  increases  at  every  point  except  the  stagnation  zone  for 
this  case.  The  on-target  irradiance  also  increases  at  every  point  except  the  stagnation  zone 
(along  the  vertical  axis  in  this  case). 


Figure  16. 


Longitudinal  Range  (km) 

Simple  case  scenario  for  =  10(m /  s)  y,  =  -3(m /  s)  y, 
hj  =  \0m,  P  =  200 kW and/l  =  1.064  pm . 
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In  Figure  17,  the  on-target  irradianee  is  plotted  for  =  10(m/ s) 

K..«.r=10(m/s)i),  y,,„,  =-3(m/s)3), /v=10m  and  P  =  200kW.  With  the 

introduction  of  target  motion  in  Figure  17,  the  cross  wind  gets  stronger  at  every  point 
except  the  stagnation  zone,  and  in  turn  the  irradianee  increases  everywhere  except  along 
this  zone. 
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Figure  17.  Simple  case  scenario  for  =  10(m  /  s)  y,  =  10(m /  s)  y, 

Kind  =  -3(m /  s)  y,  /v  =  10  m,  F  =  200 kW  and  A  =  1.064  pm . 

In  Figure  18,  the  on-target  irradianee  is  plotted  for  =10(m/s)y, 

Krget  =  55 (m /  s)  y,  Km  =  -3 (m  /  s)  y,  K  =l0m  and  P  =  200 kW.  Since  we  have  more 

target  speed  in  this  case,  we  also  have  a  stronger  cross  wind  than  other  simple  case 
scenarios.  As  a  result  of  this,  we  have  the  largest  irradianee  at  every  point  except  the 
stagnation  zone.  Since  there  is  no  change  in  cross  wind  along  the  stagnation  zone,  the 
irradianee  at  this  zone  stays  the  same  for  every  simple  case  scenario. 
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Figure  18.  Simple  case  scenario  for  =  10(m  /  s)  %  V target  =  55(m  /  s)  y. 

Kind  =  -3(m /  s) }),  /v  =  10  m,  F  =  200 kW  and  A  =  1.064  pm . 

2.  Complex  Cases 

The  on-target  irradiances  are  plotted  for  more  complex  cases  in  this  section.  In 
Figure  19,  the  on-target  irradiance  is  plotted  for  =  2(m /s).r-2(m/s)  y,  hj=l0m 

and  P  =  200  kW.  There  is  no  relative  platform  and  target  motion  for  this  plot.  Since  the 
wind  blows  from  315°  azimuth,  we  have  a  stagnation  zone  on  the  line  that  connects 
315°  and  135°  azimuth,  whereas  the  cross  wind  gets  stronger  through  the  line  that 
connects  45°  and  225° .  As  a  result  of  these  zones,  we  have  a  on-target  irradiance  pattern 
that  is  compressed  along  the  135° -3 15°  line  and  elongated  along  the  45° -225°  line. 
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Longitudinal  Range  (km) 

Figure  19.  Complex  case  scenario  for  =  2(m  /  s)  x -  2(m /  s)  /v  =  10  m, 

P  =  200  kW  and  A  =  1.064  . 

In  Figure  20,  the  on-target  irradiance  is  plotted  for  =10(m/s)5', 

=2(m/s)x-2(m/s)5',  hj=l0m  and  F  =  200kW.  With  the  introduction  of 

relative  platform  and  target  motion,  the  magnitude  of  the  effective  cross  wind  changes  at 
every  point  for  this  case.  This  results  in  a  new  orientation  of  the  on-target  irradiance 
pattern. 


Longitudinal  Range  (km) 

Figure  20.  Complex  case  scenario  for  =  10(m  /  s)  y, 

=2(m/s)x-2(m/s)y,  =  10  m  and  P  =  200kW. 
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In  Figure  21,  the  on-target  irradianee  is  plotted  for  =  10(m/ s) 

=10(m/s)x,  =2(m/s)x-2(m/s)j),  =10m  and  P  =  200kW.  With  the 

introduction  of  target  motion  in  Figure  21,  the  effective  cross  wind  changes  at  every  point 
again,  and  in  turn  this  results  in  a  new  on-target  irradianee  orientation. 


With  Platform-Taraet  Motion 
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Figure  21.  Complex  case  scenario  for  =10(m/s)5',  y,arget  =10(m/s)x, 

K.W  =2(m/s)x-2(m/s)y,  /v  =10m,  P  =  200kWand  2  =  1.064 

pm. 

In  Figure  22,  the  on-target  irradianee  is  plotted  for  1?^;^  =10(m/s)5', 
=55(m/s)x,  =2(m/s)Jc-2(m/s)y, /v  =10m  and  P  =  200kW.  With  the 

increase  in  target  speed,  both  on-target  irradianee  magnitude  and  pattern  orientation 
changes  for  every  target  position  on  the  plot. 
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Figure  22.  Complex  case  scenario  for  =10(m/s)5',  y,arget  =  55(m/s)x, 
=2(m/s)x-2(m/s)5),  /v=10m,  P  =  200 kW  and  A  =  1 .064 

pm. 

D.  PRACTICAL  EXAMPLES 

In  order  to  characterize  the  effect  of  relative  target  and  platform  motion  on  the 
atmospheric  propagation  of  HELs,  practical  examples  will  be  evaluated  via  the  changes 
in  on-target  irradiance  for  different  cases  in  this  section. 

1.  UAV  Target 

In  this  sub-section,  the  on-target  irradiances  will  be  evaluated  with  varying  UAV 
velocities,  varying  platform  (ship)  velocities,  varying  wind  velocities,  different  UAV 
altitudes  and  initial  laser  powers. 

a.  The  Effect  of  Target  Altitude  and  Target  Speed 

In  Figures  23  through  25,  the  on-target  irradiances  are  plotted  for  constant  laser 

output  power,  wind  velocity  and  platform  velocity  which  are  P  =  200  kW, 

^wind  =  -3(m  /  s)x  and  =  2(m !  s)x  +  2(m  /  S) y ,  respectively.  Target  velocities 

change  from  20(m/s)  y  to  198(m/s)  y  .  In  Figure  23,  the  on-target  irradiances  are  plotted 

for  target  height  Hj  =  300  m.  The  amount  that  the  atmospheric  propagation  is  effected 

by  the  relative  target  and  platform  motion  decreases  with  increasing  speed  of  the  target. 
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As  the  target  speed  increases,  the  effective  cross  wind  starts  to  be  dominated  by  the  target 
speed,  and  cooling  of  the  laser  beam  path  leads  to  much  less  thermal  blooming. 
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Figure  23.  On-target  irradiances  for  varying  target  velocities. 

Kind  =  -3(m  /  s) V ,  =  2(m  /  s) V  +  2(m  /  s)  y ,  P  =  200  kW  and 

=300m. 


In  Figure  24,  the  on-target  irradiances  are  plotted  for  Hj^  =  800  m.  The  same 

trends  are  seen  as  in  Figure  23.  Also,  as  the  target  altitude  increases,  the  slew  rate  of  the 
beam  director  following  the  target  decreases  and  in  turn  the  effect  of  relative  target 
motion  decreases. 
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Figure  24.  On-  target  irradiances  for  varying  target  velocities, 

y«in4  =  -3(m  /  s)x ,  =  2(m  /  s)x  +  2(m  /  s) y ,  F  =  200  kW  and 

=800m. 


In  Figure  25,  the  on-target  irradiances  are  plotted  for  =1300  m.  Again,  the 
same  trends  are  observed  as  in  the  previous  figures. 
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Figure  25.  On-target  irradiances  for  varying  target  velocities, 

y«in4  =  -3(m  /  s)x ,  =  2(m  /  s)x  +  2(m  /  s) y ,  F  =  200  kW  and 

=  1300  m. 


b.  The  Effect  of  Laser  Output  Power 

In  Figure  24,  the  on-target  irradiances  are  plotted  for  constant  target  heights,  wind 
velocities,  and  platform  velocities,  which  are  Hj  =%Q0m,  = -3(m / s)v  and 

Kft/p  =  3(m/s)v-i-2(m/s)y ,  respectively.  Target  velocities  and  laser  output  powers 

change  from  20(m/s).r  to  60(m/s)x  and  200  kW  to  1  MW,  respectively.  It  is  seen  that 
the  effect  of  the  relative  target  and  platform  motion  on  the  atmospheric  propagation 
depends  on  the  output  laser  power,  since  target/platform  motion  mainly  affects  thermal 
blooming  which  is  more  prevalent  at  higher  laser  powers. 
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Figure  26.  On-target  irradiances  for  varying  target  velocities  and  laser 

output  powers,  =  -3(m  /  s)x  ,  =  2(m  /  s)x  +  2(m/  s)}'  and 

=800m. 


c.  The  Effect  of  Target  Direction 

In  Figure  27,  the  on-target  irradiances  are  plotted  for  constant  laser  output  power, 
wind  velocity,  platform  velocity,  and  target  speed,  which  are  P  =  200kW, 

Kind  =-3(m/s)v,y^;„^  =2(m/s)v-L2(m/s)y  and  =122  m/s,  respectively.  Target 

direction  and  target  altitudes  change  from  0°  to  90°  azimuth  and  300  meters  to 
1300  meters,  respectively.  Figure  27  shows  that  the  direction  of  the  target  affects  the 
atmospheric  propagation  such  that  the  orientation  of  the  on-target  irradiance  contours 
changes  in  a  similar  manner  with  target  heading  regardless  of  the  target  height. 
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Figure  27.  On-target  irradiances  for  varying  target  directions  and  heights, 
K.W  =  -3(m  /  s)v  ,  y,  =  2(m  /  s)x  +  2(m  /  s) y ,  P  =  200  kW  and 


^ target  =  122  m/S. 


d.  The  Effect  of  Platform  Velocity 

The  effect  of  the  platform  velocity  on  the  effect  of  relative  target  and  platform 
motion  on  the  atmospheric  propagation  will  be  studied  with  different  platform  speeds  and 
platform  headings. 
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( 1 )  The  Effect  of  Platform  Direction 

In  Figure  28,  the  on-target  irradiances  are  plotted  for  constant  laser  output 
power,  wind  velocity  and  target  height,  which  are  P  =  200  kW,  =  -3(m  !  s)x  ,  and 
Hj  =500  m,  respectively.  Target  speed  changes  from  20  m/s  to  180  m/s  in  the  y- 

direction  and  platform  direction  changes  but  with  a  constant  speed  of  9  m/s.  As  the  target 
speed  increases,  it  starts  to  dominate  the  irradiance  pattern  and  the  platform  direction  has 
less  of  an  impact. 
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Figure  28.  On-target  irradiances  for  varying  target  speeds  and  platform 
headings,  =  -3(m  /  s)v ,  P  =  200  kW  and  Hj  =  500  m. 
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(2)  The  Effect  of  Platform  Speed 

In  Figure  29,  the  on-target  irradiances  are  plotted  for  constant  laser  output  power, 
wind  velocity  and  target  height,  which  are  P  =  200kW,  = -3(rr\/ s)x ,  and 

Hj  =500  m,  respectively.  Target  speed  changes  from  20  m/s  to  180  m/s  in  the  y- 

direction  and  platform  speed  changes  from  5  m/s  to  15  m/s  in  the  x-direction.  As  the 
target  speed  increases,  it  starts  to  dominate  the  irradiance  pattern  and  the  platform  speed 
has  less  of  an  impact. 
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Figure  29.  On-target  irradiances  for  varying  target  and  platform  speeds, 
'^wind  =  “3(m !  s)x ,  P  =  200  kW  and  Hj.  =  500  m. 
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e.  The  Effect  of  the  Wind  Velocity 

The  effect  of  the  wind  velocity  on  the  effect  of  relative  target  and  platform  motion 
on  the  atmospheric  propagation  will  be  studied  with  different  wind  speeds  and  wind 
headings. 

( 1 )  The  Effect  of  Wind  Direction 

In  Figure  30,  the  on-target  irradiances  are  plotted  for  constant  laser  output  power, 
platform  velocity  and  target  height,  which  are  P  =  200kW,  = -9(m/s)v,  and 

Hj  =800  m,  respectively.  Target  speed  changes  from  20  m/s  to  180  m/s  in  the  y- 
direction  and  wind  direction  changes  but  with  a  constant  speed  of  4  m/s.  As  the  target 
speed  increases,  it  starts  to  dominate  the  irradiance  pattern  and  the  wind  direction  has  less 
of  an  impact  on  the  irradiance  like  the  effect  of  platform  direction. 
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Figure  30.  On-target  irradiances  for  varying  wind  directions  and  target  speeds, 

=  -9(m /s)x,  P  =  200 kW  and  =  800  m. 


(2)  The  Effect  of  Wind  Speed 

In  Figure  31,  the  on-target  irradiances  are  plotted  for  constant  laser  output  power, 
platform  velocity  and  target  height,  which  are  P  =  200kW,  = -9(m/s)v,  and 

=800  m,  respectively.  Target  speed  changes  from  20  m/s  to  180  m/s  in  the  y- 

direction  and  wind  direction  changes  from  2  m/s  to  8  m/s.  As  the  target  speed  increases,  it 
starts  to  dominate  the  irradiance  pattern  and  the  wind  speed  has  less  of  an  impact  on  the 
irradiance  like  the  effect  of  platform  speed. 
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Figure  3 1 .  On-target  irradiances  for  varying  wind  and  target  speeds, 
=  -9(m /s)x,  P  =  200 kW  and  =  800  m. 


2.  Small  Boat  Target 

In  this  sub-section,  the  on-target  irradiances  will  be  evaluated  with  varying  target 
(small  boat)  velocities,  varying  platform  (ship)  velocities,  varying  wind  velocities, 
varying  laser  output  powers  and  constant  boat  height  {Hj  =  0.4m). 

a.  The  Effect  of  Target  Velocity 

The  effect  of  relative  target  and  platform  motion  on  the  atmospheric  propagation 
will  be  studied  with  different  target  speeds  and  target  headings. 
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(1)  The  Effect  of  Target  Direction 

In  Figure  32,  the  on-target  irradiances  are  plotted  for  constant  wind  velocity, 
platform  velocity  and  target  height,  laser  output  power,  and  target  speed,  which  are 

^w,w=-3(m/s)v,  y,,,,=5(m/s)y,  H,=QAm,  P  =  200  kW,  and  m/s, 

respectively.  The  direction  of  the  target  changes  for  each  plot.  It  is  seen  that  the 
atmospheric  propagation  is  affected  in  each  plot  with  the  change  in  target  direction. 
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Figure  32.  On-target  irradiances  for  varying  target  directions, 

k,w  =  -3(nn  /  s)i ,  V,„„  =  5(m  /  s)^ ,  ^>  =  200  kW,  =  9  m/s  and 

Hj  =  0.4  m. 


(2)  The  Effect  of  Target  Speed 

In  Figure  33,  the  on-target  irradiances  are  plotted  for  constant  laser  output  power, 
wind  velocity,  platform  velocity  and  target  height,  which  are  P  =  200  kW, 

^wind  =-3(m/s)v,  =5(m/s)5'  and  Hj  =0.4m,  respectively.  Target  speed  changes 

from  0  m/s  to  18  m/s  in  the  x-direction.  The  on-target  irradiance  plot  changes  noticeably 
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until  the  target  speed  dominates  the  relative  motion.  After  this  speed,  the  ehange  is  very 
small. 


-5  0  5 

Longitudinal  Range  (km) 


-5  0  5  -5  0  5 

Longitudinal  Range  (km)  Longitudinal  Range  (km) 


V  =5(m/s)y,  V.  =14(m/s)x 

_  S|)(P  '  "  target  '  g 


-5  0  5 

Longitudinal  Range  (km) 


V3,.=5(m/s)y.V,^  =18(m/s)x 


-5  0  5 

Longitudinal  Range  (km) 


V,,,=5(m/s)y,V,^  t=''°('"/®)" 


-5  0  5 

Longitudinal  Range  (km) 


Figure  33.  On-target  irradiances  for  varying  target  velocities, 

=  -3(m  /  S)i  ,  =  5(m  /  s)y ,  P  =  200 kW  and  H,  =  0.4  m. 


b.  The  Effect  of  Laser  Output  Power 

In  Figure  34,  the  on-target  irradiances  are  plotted  for  constant  wind  velocity,  platform 
velocity  and  target  height,  which  are  =-3(m/s)v,  =l{ml  s)x  +  l{xr\l s)y  and 
Hj  =0.4m,  respectively.  Target  speed  changes  from  0  m/s  to  18  m/s  in  the  x-direction 

and  laser  output  power  changes  from  200  kW  to  1  MW.  It  is  seen  that  the  effect  of  the 
relative  target  and  platform  motion  on  the  atmospheric  propagation  increases  with  the 
increase  in  the  output  laser  power. 
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Figure  34.  On-target  irradiances  for  varying  target  velocities  and  laser  output 
powers,  =  -3(m  /  s)x ,  =  2(m  /  s)x  +  (m/s)y  and  Hj  =  0.4  m. 


c.  The  Effect  of  Platform  Velocity 

The  effect  of  the  platform  velocity  on  the  effect  of  relative  target  and  platform 
motion  on  the  atmospheric  propagation  will  be  studied  with  different  platform  speeds  and 
platform  headings. 

(1)  The  Effect  of  Platform  Direction 

In  Figure  35,  the  on-target  irradiances  are  plotted  for  constant  laser  output  power, 
wind  velocity  and  target  height,  which  are  P  =  200  kW,  =  -3(m  /  s)x  ,  and  Hj  =  0.4 
m,  respectively.  Target  speed  changes  from  7  m/s  to  18  m/s  in  the  y-direction  and 
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platform  direction  changes  but  with  a  constant  speed  of  9  m/s.  Unlike  the  UAV  case,  the 
target  speed  does  not  dominate  the  irradiance  pattern  since  the  difference  between  target 
and  platform  speed  is  not  so  large.  As  a  result  of  this,  the  on-target  irradiance  pattern  gets 
effected  with  any  change  in  the  platform  heading,  although  this  effect  decreases  as  the 
target  speed  increases. 
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Figure  35.  On-target  irradiances  for  varying  target  speeds  and  platform 
headings,  =-3(m/s)i,  P  =  200kW  and  =0.4  m. 


(2)  The  Effect  of  Platform  Speed 

In  Figure  36,  the  on-target  irradiances  are  plotted  for  constant  laser  output  power, 
wind  velocity  and  target  height,  which  are  P  =  200  kW,  =  -3(m  /  s)x  ,  and  Hj  =  0.4 
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m,  respectively.  Target  speed  changes  from  7  m/s  to  18  m/s  in  the  y-direction  and 
platform  speed  changes  from  5  m/s  to  15  m/s.  Unlike  the  UAV  case,  the  target  speed  does 
not  dominate  the  irradiance  pattern  since  the  difference  between  target  and  platform 
speed  is  not  so  large.  As  a  result  of  this,  the  on-target  irradiance  pattern  gets  affected  with 
any  change  in  the  platform  speed,  although  this  effect  decreases  as  the  target  speed 
increases. 
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Figure  36.  On-target  irradiances  for  varying  target  and  platform  speeds, 
V^ind  =  -3(m  /  s)x ,  P  =  200  kW  and  =  0.4  m. 
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d.  The  Effect  of  the  Wind  Velocity 

The  effect  of  the  wind  velocity  on  the  atmospheric  propagation  will  be  studied 
with  different  wind  speeds  and  headings. 

( 1 )  The  Effect  of  Wind  Direction 

In  Figure  37,  the  on-target  irradiances  are  plotted  for  constant  laser  output  power, 
platform  velocity  and  target  height,  which  are  P  =  200  kW,  =  -3(m  /  s)x ,  and 
Hj  =0.4  m,  respectively.  Target  speed  changes  from  7  m/s  to  18  m/s  in  the  y-direction 

and  wind  direction  changes  but  with  a  constant  speed  of  7  m/s.  Unlike  the  UAV  case,  the 
target  speed  does  not  dominate  the  irradiance  pattern  since  the  difference  between  target 
and  wind  speed  is  not  so  large.  As  a  result  of  this,  the  on-target  irradiance  pattern  gets 
affected  with  any  change  in  the  wind  heading,  although  this  effect  decreases  as  the  target 
speed  increases. 
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Figure  37.  On-target  irradianees  for  varying  target  speeds  and  wind  headings, 

=  -3(m I s)x,  P  =  200 kW  and  =  0.4  m. 


(2)  The  Effeet  of  Wind  Speed 

In  Figure  38,  the  on-target  irradianees  are  plotted  for  constant  laser  output  power, 
platform  velocity  and  target  height,  which  are  P  =  200kW,  =-3(m/s).v,  and 

=0.4  m,  respectively.  Target  speed  changes  from  7  m/s  to  18  m/s  in  the  y-direction 

and  wind  speed  changes  from  0  m/s  to  10  m/s.  Unlike  the  UAV  case,  the  target  speed 
does  not  dominate  the  irradiance  pattern  since  the  difference  between  target  and  wind 
speed  is  not  so  large.  As  a  result  of  this,  the  on-target  irradiance  pattern  gets  affected  with 
any  change  in  the  wind  speed,  although  this  effect  decreases  as  the  target  speed  increases. 
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Figure  38.  On-target  irradiances  for  varying  target  and  wind  speeds, 
=  -3(m /s)x,  P  =  200 kW  and  =  0.4  m. 
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VI.  CONCLUSIONS 


According  to  the  results,  the  effects  of  relative  target  and  platform  motion  on  the 
atmospheric  propagation  of  HELs  depend  on  parameters  like  target  altitude,  target 
velocity,  platform  velocity,  laser  output  power,  and  wind  velocity.  As  the  target  altitude 
increases,  the  effect  of  the  relative  motion  on  the  propagation  decreases,  because  the  slew 
rate  of  the  beam  director  decreases  with  the  increase  in  target  altitude.  The  wind  velocity 
can  either  decrease  or  increase  the  relative  motion  effect  on  propagation  for  a  specific 
target  acquisition,  depending  on  the  specific  velocities  of  the  platform  and  target  as  well 
as  their  relative  positions.  As  an  example,  as  the  target  speed  becomes  faster  than  the 
platform  speed,  the  relative  motion  becomes  dominated  by  the  target  speed.  Changes  in 
both  target  and  platform  direction  affect  the  propagation  of  the  laser  regardless  of  the 
target  altitude,  target  and  platform  speeds,  and  the  wind  velocity,  because  the  effective 
cross  wind  changes  with  the  change  in  directions  for  each  case.  The  effect  of  the  relative 
target  and  platform  motion  on  the  atmospheric  propagation  also  depends  on  the  output 
laser  power,  since  target/platform  motion  mainly  affects  thermal  blooming,  which  is 
more  prevalent  at  higher  output  laser  powers. 

As  a  result,  the  relative  target  and  platform  motion  significantly  affects  the 
atmospheric  propagation  of  HELs,  and  the  level  that  the  propagation  is  affected  depends 
on  the  parameters  discussed  above.  HELs  can  be  used  more  effectively  by  the  users 
considering  these  effects. 
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